Fused Azulenes: Possible Organic Multiferroics 



Simil Thomaa3 and S. RamaseshgO] 

Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore 560 012, India 

Karen Hallber^J and Daniel Garcic^j 
Centra Atomico Bariloche and Instituto Balseiro, 
Comision Nacional de Energia Atomica and CONICET, 8400 Bariloche, Argentina 

(Dated: November 30, 2011) 

We study low-lying correlated electronic states of fused azulenes using the long-range interacting 
Pariser-Parr-Pople model and the finite density matrix renormalization group method. The ground 
state is a singlet for oligomers with up to 5 azulene units. For oligomers with more than 5 azulene 
units and beyond 6 azulene units the oligomers have a triplet ground state. From the excitation gaps 
between the lowest Ms = state and the lowest states in Ms = 1, 2 and 3 sector we predict that the 
ground state spin of the fused azulene increases with number of azulene units. In the thermodynamic 
limit we expect the fused azulene to be a ferromagnet. Charge density calculations show that the 
ground state of the system has ferroelectric alignment of the dipoles of the monomeric units. Thus, 
a fused azulene system could be the first example of a molecule which is both ferromagnetic and 
ferroelectric in the ground state. 

PACS numbers: 71.10.Fd, 75.50.Xx, 77.55.Nv, 77.84.Jd 



INTRODUCTION 

Conjugated organic materials have been attracting 
much attention for many electronic applications such as 
light emitting diodes[U [2], photovoltaic cells(PVC)[3], 
and thin film transistors [1], to name a few. These ma- 
terials are usually conjugated polymers such as MEH- 
PPV0, PF2/60, PHT[7] which comprise of conjugated 
back bones built from phenyl or other ring systems. In 
recent years, fused systems such as oligo acenes have also 
come to the fore [8HTT]. Pentacene has become impor- 
tant because its lowest triplet state has an excitation 
energy which is slightly less than half the singlet exci- 
tation gap, allowing multi exciton generation (MEG)[TT] 
through singlet fission. If the open cell voltage in triplet 
dissociation is not small compared to singlet dissociation, 
the efficiency of the PVC can be significantly enhanced. 
This new paradigm has inspired a search for conjugated 
systems with low energy triplets. 

One design principle that can lead to a low energy 
triplet state is incorporation of frustration in transfer 
topology of the molecule. The ground state as well as 
the lowest triplet state are dominated by covalent states 
(states with one electron per site). Their stabilization 
arises from derealization of electrons in these states. In 
the singlet state the derealization is reduced if the sys- 
tem is frustrated thereby increasing its energy. However, 
the triplet state in frustrated systems does not experience 
such enhanced frustration. This results in a reduction of 
the singlet-triplet gap. 

A simple frustrated conjugated molecule that has been 
well studied is azulene with fused five and seven member 
rings. It also has a dipole moment ~ l.OD attributed 
to the Hiickel 4n + 2 rule[T2]. Fused azulenes have long- 



range frustration and could lead to low-lying triplets. Be- 
sides, the dipole moment of the azulene monomer can also 
align resulting in a molecular ferroelectric phase. With 
this in focus, we have carried out a detailed study of the 
low-lying excitations of fused azulene using finite den- 
sity matrix renormalization group (DMRG) technique. 
In the next section, we briefly discuss the model Hamil- 
tonian and the DMRG method used in our study. In the 
last section we present our results and discussions. 

MODEL HAMILTONIAN AND THE DMRG 
METHOD FOR FUSED AZULENES 

The 7r system in fused polyazulenes (Fig. [T^,) is mod- 
eled by considering one p z orbital on each carbon atom. 
The Pariser-Parr-Pople (PPP) model pi Hi] Hamilto- 
nian is given by, 

Hppp = ^2 *Ol,<A a + H - c -) + -^^niifii - 1) 
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with the symbols having their usual meanings. The pen- 
tagons and the heptagons are treated as regular polygons 
of side 1.397 A. The transfer integral between bonded 
sites is taken to be -2.4 eV. The Hubbard parameter U 
for Carbon is fixed at 11.26 eV. The inter-site interactions 
are parameterized using the Ohno[T5] formula. The PPP 
model with these standard parameters have been exten- 
sively used for successfully modeling the excitations in a 
host of conjugated molecules. 

The DMRG [16] method initially developed was a tech- 
nique for studying correlated Hamiltonians such as the 
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Hubbard and Heisenberg models. In chemistry, the 
method has proved to be ideally suited to study conju- 
gated polymeric system, within model Hamiltonian ap- 
proaches. It has been shown that the DMRG method is 
accurate even in the presence of long-range interactions 
if these interactions are diagonal in real space, as with 
PPP model. This follows from comparing the DMRG re- 
sults with exact diagonalization results for short chains 
as well as from entanglement entropy analysis of their ex- 
act eigenstates[T7]. Many of the interesting conjugated 
polymers consist of ring systems, such as the phenyl, pyri- 
dine, thiophene, furan or pyrrole rings. Building the de- 
sired oligomer by adding two sites at each infinite DMRG 
step is nontrivial and the general method was first illus- 
trated for the case of poly para phcnylcncs 18J. A similar 
approach can be used for building fused polyazulenes. 
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FIG. 1. (a)Structure of fused azulene, (b)DMRG scheme for 
building fused azulene. 

The starting point in the infinite DMRG scheme is a 
four membered ring, at each iteration two sites are added 
in the middle of the ring and the ring size is expanded till 
we reach a ring size of ten at which point, an additional 
bond is introduced between the new site and a site added 
at the previous iteration, both being in the same block 
(Fig. [T] b ) . The procedure is continued until we reach 
the oligomer of the desired size. After this point we start 
the finite DMRG sweeps in the usual manner [16|. 

We have checked our DMRG calculations against the 
noninteracting model results for all oligomers, for differ- 
ent cut-offs in the number of dominant density matrix 
eigenvectors, m. The ground state energy for a system of 

3 fused azulenes agrees to almost Ippm for m = 500. We 
have also checked the accuracy of our results for smaller 
system sizes and have found that the the ground and 
as the lowest excited state energies agree to better than 

4 x 10~ 4 eV for this cut-off. Hence in all our studies we 
have set m — 500. 



RESULTS AND DISCUSSION 

Fused azulenes can have different structures, in one 
structure, the upper "polyene" chain is trans-trans while 
the lower "polyene" chain is cis-cis(Fig. [l^,). In another 
structure, both the upper and lower "polyene" chains are 
cis-trans. We have carried out the DMRG calculations on 
the structure in Fig. [T^,. Since the difference between the 
two structures arises only in third nearest neighbor, the 
electronic structure is not expected to be very different 
in the two cases. The energy per azulene unit linearly 
extrapolate to —19.117 eV in the thermodynamic limit. 

Spin gap in these oligomers is usually computed as the 
difference in energy between the lowest energy states in 
the M s = 1 sector and that in the M s = sectors (the 
system has an even number of electrons). This is based 
on the assumption that the ground state is a singlet and 
the lowest spin excited state is a high spin state; the sin- 
glet state is found in the M s = sector and states with 
nonzero spin are found in the M s = 1 sector. Further- 
more, it is often assumed that the lowest energy spin state 
in the M s = 1 space is a triplet. Accordingly, we com- 
puted the lowest energy states in the M s = and M s = 1 
subspaces. Our studies showed that for oligomers with 
more than five azulene units the lowest energy states in 
these two subspaces were degenerate within the accuracy 
of the calculation. This implies that the ground state has 
a spin S > 1. 

To determine the spin of the ground state in the longer 
oligomers (n > 6), we computed two lowest energy states 
in the M s = and M s = 1 subspaces and the lowest en- 
ergy state in the M s = 2, 3 subspace. In these oligomers, 
the lowest energy state in M s = 2 state is higher in en- 
ergy than the degenerate lowest energy states in the M s 
= and 1 subspaces. This conclusively proves that the 
ground state in these oligomers is a triplet. To investi- 
gate whether in even longer fused azulenes, the ground 
state again switches from S = 1 to higher spin, we in- 
vestigated the gaps between the lowest energy Ms = 
state and lowest energy states in Ms = 1,2 and 3, as 
function of the number of oligomers (Fig. [2]). We note 
that the slopes for these different gaps in Fig. [2] are dif- 
ferent and it appears that as the system size increases, 
the ground state will switch initially from S = to S — 1 
as seen from the vanishing of the gap between the lowest 
energy levels in the Ms = 1 and Ms = sector. From 
the behavior of the gap between the lowest Ms = 2 and 
Mg = level, we see that these gap is also going to vanish 
when the system size increases to about 10 or 11 azulene 
units leading to spin S = 2 ground state. The gap A 3 
between lowest Ms — 3 and Ms — levels also appears 
to vanish for even larger systems, making 5 = 3 ground 
state. Therefore, it is likely that in the polymer limit, 
a ferromagnetic ground state will result. Polymers with 
triplet ground state are, to the best of our knowledge not 
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FIG. 2. Spin gaps (in eV) vs ±. ( Ai = (£W s =i - £a/ s =o ), 
A 2 = (E Ms =2-E Ms =o), and A 3 = (£m s=3 -B« s =o). These 
states correspond to the lowest energy state in their respective 
M a sector. 



FIG. 3. Ring charges in fused azulenes with 8, 9 and 10 azu- 
lene units. Charge in a ring is given by (n — p), p =number of 
carbon atoms in the ring. Odd(Even) ring index corresponds 
to seven(five) membered ring of the azulene. 



known, although Monkman and coworkers |19j. based on 
the correlation between experimental singlet and triplet 
excitation energies, predicted that a polymer will have a 
triplet ground state, if the singlet exciton energy is be- 
low 1.3 eV. However, there is no theoretical basis for this 
prediction and is just based on extrapolation of the plot 
of singlet exciton energy vs triplet exciton energy from 
experiments, for various polymers. To see if the system 
is a polymeric example of flat-band ferromagnetism, we 
examined the non-interacting one-particle states of the 
system with 12 azulene units. We find that there is no de- 
generacy of the nonbonding molecular orbitals(MOs) at 
Ep, which is a requisite for flat band ferromagnetism [2D]. 
These results show that for fused azulene with less than 
six azulene units, the MEG process can create multiple 
triplet excitations. However, it is not clear whether it 
will lead to higher efficiencies in a PVC as the voltage 
obtained for the cell could be quite low, due to small 
triplet gaps. 

To understand whether the high spin ground state re- 
sults from isolated azulene units in the polymer or is a 
specific property associated with fused systems, we have 
also carried out DMRG studies on polyazulene which has 
two azulenes per unit cell and has successive azulene 
units connected by a C — C bond between pentagons. 
We have studied the spin gap in these systems by us- 
ing finite DMRG algorithm with a cut-off in m of 400. 
For all sizes we have studied, we find that singlet is the 
ground state. The energy gap to the triplet state lies 
in the range 1.08 ± 0.04 eV for the oligomers with 4, 6 
and 8 azulene units. Thus, for observing a high spin 
ground state, it appears that the azulenes have to be 



fused so that the bond frustration extends over the en- 
tire system. The polyazulene system has only localized 
bond-frustration, hence the singlet ground state is pre- 
ferred. Bond-order of a bond (i,j) in the ground state is 
defined as \ < ^2 a {al (7 a j.cr+ a j a a i,cr} > and expectation 
value is the ground state expectation value. Bonds with 
large (small) bond orders will have a tendency to become 
shorter(longer) in the equilibrium geometry and thus in- 
dicate the equilibrium bond lengths when the system is 
allowed to relax its geometry. We note that in both the 
singlet and the triplet states, the bond orders for the pe- 
ripheral bonds are nearly uniform. This shows that the 
system is unlikely to undergo a structural distortion. 

This observation is important, since it is well known 
that low-dimensional systems can undergo a distortion 
of the Peierls' type which could lead to stabilization of 
the singlet state relative to the triplet state. We also 
note that the bond orders for the bonds that connect 
the upper and lower polyene like chains is rather small, 
implying that any distortion of the system will only uni- 
formly increase the distance between the upper and lower 
chains, which is unlikely to alter energy level ordering of 
the low-lying states. 

We have investigated the charge distribution in the 
oligomers. We have shown the net charge (equals sum of 
actual charges on the the sites comprising the ring minus 
the number of sites in the ring) in the seven and five mem- 
bered rings in the oligomers in Fig. [3| In all the oli gomcrs 
studied, seven membered rings have positive charge while 
five membered ring have negative charge. The amplitude 
of the charge density wave in the system remains approx- 
imately the same, independent of the oligomer size and 
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FIG. 4. Charge gap (in eV) vs 



we expect a ferroelectric state in the polymer limit. This 
would mean that the fused azulene system in the ther- 
modynamic limit, would be a multiferroic, showing both 
a ferromagnetic and a ferroelectric ground state. 

Charge gap measures the energy required to create an 
independent electron-hole pair in a system. The charge 
gap for a polymeric system can be obtained by extrapo- 
lating E C (N) = E+(N) + E~(N) - 2E° S (N) (where 
N is the number of monomers and superscripts '+' and 
'— ' refer to cation and anion while '0' refers to the neu- 
tral species) to the thermodynamic limit. The quantity 
E C (N) gives the energy required to create a pair of free 
moving electron and hole in the ground state. The plot 
of charge gap vs ^ is shown in Fig. |4j Extrapolating the 
charge gap to thermodynamic limit, we obtain the value 
Ec{oo) = 0.403 eV. This is smaller than typical charge 
gaps found for conjugated polymers by almost an order 
of magnitude pi]. 

To conclude, we have studied low-lying electronic 
states in fused azulenes using PPP model and finite 
DMRG scheme. We note that for short oligomers (< 6 
azulene units) the ground state is a singlet and the 
ground state switches to triplet for all the fused systems 
between 6 and 10 azulene units long. In the thermody- 
namic limit, we expect the chain to be ferromagnetic. 
This implies that MEG process could be important only 
in short fused azulene systems where the triplet is a low 
energy excited state. Bond order studies show that the 
bond along the upper and lower chain will remain uni- 
form while it is likely that the separation between the 
upper and lower chains of fused azulene will increase in 
thermodynamic limit. The system also exhibits order- 



ing of dipoles of the azulenes with seven membered rings 
being positively charged and five membered rings nega- 
tively charged. This implies that the ground state will 
be both a ferromagnet and ferroelectric in the thermody- 
namic limit. Such a system would be the first example of 
a molecular multiferroics. The charge gap in these sys- 
tem is quite small compared to other polymers. These 
features of fused azulene should be of importance in or- 
ganic device applications. 
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